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ABSTRACT 

With the temperature power spectrum of the cosmic microwave background (CMB) at least 
four orders of magnitude larger than the 5-mode polarisation power spectrum, any instrumen¬ 
tal imperfections that couple temperature to polarisation must be carefully controlled and/or 
removed. Here we present two new map-making algorithms that can create polarisation maps 
that are clean of temperature-to-polarisation leakage systematics due to differential gain and 
pointing between a detector pair. Where a half wave plate is used, we show that the spin-2 
systematic due to differential ellipticity can also be removed using our algorithms. The al¬ 
gorithms require no prior knowledge of the imperfections or temperature sky to remove the 
temperature leakage. Instead, they calculate the systematic and polarisation maps in one step 
directly from the time ordered data (TOD). The first algorithm is designed to work with scan 
strategies that have a good range of crossing angles for each map pixel and the second for scan 
strategies that have a limited range of crossing angles. The first algorithm can also be used to 
identify if systematic errors that have a particular spin are present in a TOD. We demonstrate 
the use of both algorithms and the ability to identify systematics with simulations of TOD 
with realistic scan strategies and instrumental noise. 

Key words: methods: data analysis - methods: statistical - cosmology: cosmic microwave 
background - cosmology: large-scale structure of Universe 


1 INTRODUCTION 

The CMB contains an incredible wealth of cosmological informa¬ 
tion. The properties of the Universe can be probed at a number 
of different epochs using the CMB. The very early universe can 
be probed through the CMB’s constraints on inflation parameters 
( [Planck Collaboration et al.|2014d| l. Physics before last scattering 
is imprinted on the CMB as baryon acoustic oscillations, these os¬ 
cillations have been mapped to exquisite detail using both the tem¬ 


perature power spectru m ([Planck Collaboration et al.|2014b[[Stor^ 


[et al.[[2013|[Das et al.[ 

[2014 

and the 5-mode polarisation power 

spectrum ([Crites et al. 

|2014 

[Brown et al.[[2009b| [Planck Collab-] 


[oration et al.|2015a| l. The large-scale structure of the universe can 

also be probed via gravitational lensing of the CMB. This effect 
has been measured using high resolution temperature maps of the 


CMB ( [Baxter et al. [[20141 [Das et ar|[2014[ [Planck Collaboration] 

[et al.|201^. 

The CMB 5-mode polarisation power spectrum contains ad¬ 
ditional information on two of these epochs. 5-mode polarisation 
on large angular scales provides us with the best insight into in¬ 
flation by placing a direct constraint on the tensor-to-scalar ratio. 
Tentative measurements have been made in this area by BICEP2 
( [Ade et al.|2014[ ). However, Galactic foreground emission from po¬ 
larised dust has been shown to be responsible for some and pos- 
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sibly all of the signal detected (Planck Collaboration et al.|2014a 


[BICEP2/Keck and Planck Collaborations et al.[[2015[ l. The small 

scale 5-mode power spectrum is a result of gravitational lensing of 
the larger 5-mode power spectrum, which has in recent years been 
detected by a number of experiments (The Polarbear Collaboration:] 


[P. A. R, Ade et al. [2014[ [Hanson et al.|2013[ l 


As the 5-mode power spectrum is at least four orders of mag¬ 
nitude below the temperature power spectrum any coupling be¬ 
tween the two signals due to instrumental imperfections must be 
carefully controlled and/or removed. Approaches used in the lit¬ 
erature to ensure the validity of a polarisation map broadly fall 
within two categories. The first of these relies on detailed simu¬ 
lations of the instrumental set up and uses knowledge of the tem¬ 
perature sky to simulate the effects of any imperfections on the re¬ 
covered 5-mode power spectrum. This was done very successfully 
by the POLARBEAR collaboration in their detection of the lensing 
5-modes ( [The Polarbear Collaboration: P A. R. Ade et sl \.\2014\ . 
The second category involves calculating the coupling by fitting the 
parameters describing the imperfections, and then subtracting this 
coupling using CMB temperature measurements. This technique 
was shown to be effective in the analysis of BICEP2 ( [Ade et ^ 
[2014] see their figure 5). However, there is a question as to whether 
this de-projection technique would work as effectively with a more 
complex scan strategy. In addition, the fitting procedure employed 
also removes some polarisation signal. This results in a leakage 
of 5-modes to 5-modes which must be simulated and removed in 
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the power spectrum estimation ( |BICEP2 Collaboration et al.|2015^ . 
Here we present alternative novel algorithms to identify and remove 
some of the systematics that are problematic in CMB polarisation 
experiments. A key feature of our approach is that it does not re¬ 
quire any prior knowledge of the telescope or CMB temperature 
field. In addition, since there is no fitting involved, our techniques 
do not result in any leakage of ^-modes to B-modes. 

[Wallis et^ ( |2014j ) suggested a map-making algorithm to re¬ 
move systematics for experiments where there is no half-wave- 
plate (HWP). The method consists of two stages; first systematics 
of a different spin to those we want to measure are removed (spin-0 
for temperature and spin-2 for polarisation), then a second cleaning 
procedure is required to remove systematics of the same spin. |Wal-1 
jlis et al.| ( |2014| ) concentrated on beam systematics. Consequently, 
the potential source of spin-2 systematics that could couple temper¬ 
ature to polarisation that they considered was the second azimuthal 
mode of the temperature beam. To remove this they required knowl¬ 
edge of the beam to correctly predict this leakage from a tempera¬ 
ture map. This method is similar to that used by jAde et al.| ( |2014'] ) 
to remove temperature to polarisation leakage from differential el- 
lipticity. 

One potential complication that the approaches just described 
suffer from is the requirement to correctly characterise the elliptic- 
ity of the beam and hence predict the resulting leakage to polarisa¬ 
tion. Conversely, the novel approaches that we present here require 
minimal knowledge of the nature of the leakage. The methods are 
appropriate for differencing experiments that use a stepped or rotat¬ 
ing HWP. One of the algorithms can be used even if a HWP is not 
present. However, in this case only systematics of a different spin 
to polarisation can be removed. Where a HWP is present, it can be 
used to disentangle the spin-2 leakage from temperature to polar¬ 
isation due the ellipticity of the beam and the spin-2 polarisation 
signal. The two methods differ in the scan strategies to which they 
can be applied. One algorithm is suited to scan strategies where 
each map pixel is seen at a range of telescope orientations, for ex¬ 
ample the proposed EPIC scan strategy ( [Bock et ^|2009| l. The 
other is suitable for experiments where the range of orientation an¬ 
gles for each pixel is limited, for example the LSPE scan strategy 
( jAiola et al.|2012j ). There is no reason why the two methods can¬ 
not be used on different portions of the same map. Eor example 
the Planck scan strategy ( [Planck Collaboration et al.|2014b| ) results 
in good orientation coverage at the ecliptic poles where the first 
method would be most suited and a limited range at the ecliptic 
plane where the second method would be more appropriate. 

We demonstrate that our techniques can also be used to re¬ 
move differential gain and pointing even in the absence of a HWP 
with a suitable scan strategy. This does leave coupling caused 
by differential ellipticity as without a HWP this coupling is irre¬ 
ducible. In this case we advocate the previous methods of [Wallis j 
[et al.[ ( [20T4] ) and [Ade et al.[ ( [2014'] ) to remove this leakage. 

Our techniques for removing systematics involve using a 
model for the spin of the systematics, which is employed during 
map-making. This provides us with Q and U maps that are free of 
the systematics included in the model, but also maps of the sys¬ 
tematics themselves. We also demonstrate that our approach can 
be a useful method for identifying if a systematic is present in an 
experiment, or not. 

The paper is organised as follow. Section [^describes the an¬ 
alytical framework for the algorithms to remove systematics. Then 
in Section we demonstrate the use of the algorithms on simula¬ 
tions, using realistic scan strategies and time ordered data (TOD) 
which include instrumental noise. Section|^explains how the algo¬ 


rithm can be used to find systematics and demonstrates the tech¬ 
nique on a simulated TOD. Einally in Section|^we summarise our 
work. 


2 MAP-MAKING ALGORITHMS 

Our objective is to create maps free of systematic error due to the 
imperfections in the instrumentation of an experiment. We assume 
the HWP is ideal and situated at the end of the optical system (in 
emission). The effect of the HWP is to simply rotate the angle 
of the polarisation sensitivity of the beam (see, e.g. [Brown et al.[ 
[2009a| l, leaving the shape of the polarisation intensity and temper¬ 
ature beams unchanged. The assumption of an ideal HWP is obvi¬ 
ously not entirely realistic. However, in practice a HWP would only 
ever be included in an experiment if the systematic effects that they 
introduce are smaller than the effects that they are designed to mit¬ 
igate. Relaxing the assumption of HWP ideality is something we 
leave to further work. We consider an experiment where a detector 
pair is used to measure temperature and polarisation. Each detector 
is sensitive to orthogonal polarisation directions. The two signals 
di and d 2 are summed and differenced: 


= ^(di+d2), 

(1) 

= \{d,-d2). 

(2) 


In an ideal experiment would correspond to the temperature 
of the pixel and the rotated polarisation, the only effect of the 
beam would be to isotropically smooth the temperature and polari¬ 
sation fields. 

We will concentrate on recovering the polarisation of the pixel 
and therefore, we drop the superscript in equation at this stage. 
Therefore, the differenced signal, S, will be the rotated polarisa¬ 
tion of the pixel plus any systematics, the most serious of which 
will couple temperature to polarisation. Some common systemat¬ 
ics include differential gain, differential pointing and differential 
ellipticity of the detector pair. These systematics transform as spin- 
0, spin-1 and spin-2 respectively with telescope orientation - For 
a demonstration of the leakage angular dependence see e.g. fig 2 of 
[Shimon et al.[p00^ where the authors depict the monopole, dipole 
and quadrupole nature of the different systematics. The differential 
gain is spin-0 as it is simply a scaled temperature map. The differ¬ 
ential pointing is spin-1 as the signal is a difference of the temper¬ 
ature map at two close points in space. The leakage, to first order, 
is therefore proportional to the derivative of the beam smoothed 
temperature field. Differencing two elliptical Gaussians results in a 
quadrupole pattern. This quadrupole pattern is then convolved with 
the sky to create a spin-2 systematic effect. The systematic errors 
must be constant for this map-making algorithm to be able to re¬ 
move them. If the systematics change with time, a more adaptive 
algorithm would need to be developed. 

With the HWP at the end of the optical system there is no 
dependence of these systematics on the orientation of the HWP 'iph- 
The detected differenced signal S is therefore. 


5(V>h,V’t) 


,(3) 

5(V>r,V’t) 

= 3? +G + , 

(4) 


where P is the complex representation of the polarisation of the 
pixel in terms of the Stokes parameters, P = Q + ill. G, M and 
E are the temperature to polarisation leakage due to differential 
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gain, differential pointing and differential ellipticity respectively, 
and 5R is the real part operator. The magnitudes and phases of the 
systematics are dependent on the nature of the imperfections and 
the underlying temperature field. Note however that the magnitudes 
and phases are unimportant for this work. Here, only knowledge of 
the way they transform with the telescope orientation is required 
in order to remove the systematics. In equation ^ we have made 
a coordinate transformation ^r=^t+2^h. We do this so that the 
polarisation and systematics are dependent on different variables in 
our space. 

The aim of this work is, therefore, to obtain an unbiased es¬ 
timate of P, given that the detected signal depends on the system¬ 
atics as well as polarisation. The two techniques which we present 
differ only in the scan strategies that they can be applied to. We 
first present an algorithm suitable for a scan strategy where the 
'ipt coverage of a pixel is extensive. For example, the EPIC ( |Bock| 
|et al.|2009| l strategy is designed to maximise this coverage. We then 
present a second method where the coverage is limited. Balloon 
borne experiments such as LSPE ( |Aiola et al.|2012| | will have lim¬ 
ited coverage. Such experiments often include a rotatable HWP 
in order to obtain multiple polarisation crossing angles. 

When a rotating or stepped HWP is included in an experiment, 
whatever the scan strategy, the HWP can be used to provide enough 
polarisation angle coverage such that detector differencing is not 
required. As the differencing seems to lead to the temperature to 
polarisation leakage considered in this work, one may ask if other 
techniques, which do not require differencing, could be used. If the 
HWP is continually rotating then certain “lock-in” techniques can 
be used to isolate the polarisation signal from the systematic errors 
jWu et al.|2007) . However, maintaining continuous rotation of the 
HWP can cause its own wealth of systematic errors. We therefore 
focus of the case of a stepped HWP for which “lock-in” techniques 
are not applicable. 

Even with a stepped HWP, the large amount of polarisation an¬ 
gles provided by the HWP in principle allows one to recover maps 
of the Stokes parameters from just one detector. Such a technique is 
however more problematic than differencing as the temperature to 
polarisation leakage could potentially be much worse. A differenc¬ 
ing experiment allows two detectors, that are located at exactly the 
same position in the focal plane (and therefore observe the same 
point on the sky) to be used to directly remove the temperature sig¬ 
nal (see equation!^. If a single detector was used to reconstruct 
Stokes parameter maps, then the absolute pointing error, which 
is typically larger than the differential pointing considered here, 
would create different temperature responses between different ob¬ 
servations of a pixel and this would leak temperature fiuctuations 
to polarisation. 

A similar argument holds for ellipticity; by differencing de¬ 
tector pairs, we are susceptible to the difference in the ellipticity of 
two beams which often have very similar beam shapes. By creating 
polarisation maps from one detector the total ellipticity would cre¬ 
ate different temperature responses when the telescope observes a 
pixel at different orientations, leading to a much larger temperature 
to polarisation leakage. One problem that detector differencing can 
suffer from, and that using one detector avoids, is a constant differ¬ 
ential calibration. However, in this case the benefits of differencing 
often outweigh this particular problem. Another benefit in differ¬ 
encing two detectors, is that correlated noise between the detectors 
is removed. Motivated by these considerations we have adopted a 
map-making scheme that differences two detectors in a detector 
pair. 


2.1 Map-making algorithm with extensive 'ipt coverage 


Our experimental model consists of sampling a pixel at a wide 
range of orientations of the telescope and HWP. The detected sig¬ 
nal, S^, can be expressed as. 


S\pJ„ipt) 

= /i('0r,'0t)5'('0r,'0t), where 

(5) 


-j -^hits 


h{A,'ipt) 

= iVi^. E ipt) 

^ ® 7 = 1 

(6) 


is the window function representing the knowledge that we have 
of the pixel. One sample, z, will contribute one delta function 
(5('0r — 'ipl)S{'ipt — 'ipt) to this window. Our aim is to obtain an 
unbiased estimate of the polarisation of the pixel given that the sys¬ 
tematics are present and have the functional form outlined in equa¬ 
tion 0. This functional form lends itself to be described well by 
a Eourier series. We replace each term in equation with their 
Eourier series such that. 


E 




i{ni'4)r+rnx'ft) 


^ ^ hn2,m2^ 


i(n2'0r+m2'0t) 


,777-3 ^ 


iins'fr+rns'ft) 


77-2^772 

7737773 


(7) 


Multiplying each side by ^ integrating over the 

whole (^r, ^t) space and evaluating the resulting Kronecker delta 
function, we find 


Qd 

*^77 1 7771 


1 

Stt^ 


^ ^ ^772,7772*^773,7773 
7727772 
nsms 


pAtt p 27 T 

X dip, diPt e*t(’^2+n3-ui),^, + (m2+m3-mi).^t)]^(g) 
Jo Jo 


- ^ ^ ^771 —773,7771 —7773 *5^77 


nsms 


(9) 


In principle, an unbiased estimator for the different compo¬ 
nents of the signal can now be formed by inverting equation 0. 
However, this operation is not yet possible for two reasons. Eirstly, 
we are attempting to invert a matrix infinite in size. Secondly for 
any realistic window functioij^the matrix will be singular. By un¬ 
derstanding the dependence on and 'ipt of 5'(^r,'0t), we can 
ignore terms in equation (|^ where A§n 3 ,m 3 = 0 , thereby making the 
operation invertible and obtaining an unbiased estimate of Sns,m 3 
from our detected If we know the differenced signal con¬ 

tains temperature to polarisation leakage from differential gain we 
would include the term So, 0 ’ For differential pointing and ellip¬ 
ticity, we include So,±i and a§o ,±2 respectively. In principle we 
could remove systematics of any spin by simply including the cor¬ 
rect term. The polarisation of the pixel will be, 

Q = 2U{S2,o), (10) 

U = 2^(52,o). (11) 

Equation (|^ will only be invertible if there are enough hits on the 
pixel at a sufficient variety of crossing angles and HWP angles 
'iph- The more terms we include in equation the more observed 
orientations will be required. 


^ By realistic we specifically mean any window function where 
= 1. 
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2.2 Map-making algorithm with limited coverage 

The second class of experiments that we consider has a limited 
range of crossing angles and obtains polarisation angle coverage 
using a stepped HWR This is similar to the observation strategy en¬ 
visaged for the LSPE ( |Aiola et al.|2012| l. In this case using Fourier 
terms to describe the systematics is not a good choice. Here we 
describe a formalism that is specifically designed for a small, but 
non-zero, range of crossing angles. 

We start from the same position as for the case of extensive 
'ijjt coverage in Section |2.1| In an experiment we have a function 
describing the detected signal given by equation However the 
range of angles is small. This restricted range of angles means 
that describing the full Fourier mode of each systematic would be 
problematic. Instead, we choose to describe the summed effect of 
the systematics in terms of Legendre polynomials. Let the an¬ 
gles range from to We can now define a coordinate that 
spans this range: 

( 12 ) 

^max _ ^min 

where x ranges from —1 to 1. With this definition we can now 
rewrite equation 0 as 

S{^r,x) = 3 ? +/(*), ( 13 ) 

where f{x) is a function that describes the combined effects of the 
systematic leakage from temperature to polarisation. If the range 
is small enough then f(x) will be well described by only a few 
Legendre polynomials. In Fig. we show the effectiveness of the 
Legendre polynomials to describe a particular section of a function 
of the form g{'ip) = cos(2'0 + tt/S), where the '0 range is ±0.25 
rads. This range is chosen to approximate the range of crossing 
angles seen in typical balloon experiments. In particular, the maxi¬ 
mum 'ipt range in any pixel in the LSPE scan strategy ( |Aiola et ^ 
|2012| ) is ?^0.5 rads. The left panel shows that the amplitudes reduce 
almost exponentially with the order of the polynomial. In the centre 
and right panels we demonstrate that using only the first 3 Legendre 
polynomials, we can reconstruct the systematic to within fractions 
of a percent. 

With this motivation we follow similar steps to those in Sec¬ 
tion [TT] to create an estimate of the polarisation of a pixel free of 
systematics. We use the Legendre polynomials to describe the x 
dependence of the signal and a Fourier series to describe the 
dependence. To begin we write the problem as a multiple of the 
underlying signal and the window function, 

S'^{2pr,x) = h{'lpr,x)S{2pr,x). (14) 

As before we substitute the functions for their decompositions into 
a set of basis functions, where here we have chosen the Legendre 
polynomials: 

n\m\ 

Y, ( 15 ) 

n2m2 

nsms 

Taking the scalar producQof both sides with a basis function leaves 


^ The scalar product we use is, 


pAtt rl 

J d'lpr J dxf{'ipt,x)g{'ip, 


( 16 ) 


us with the triple integral, 

nd 2mi ± 1 


Stt 


^ ^ hn2,m2^n3,n 


n2m2 

nsms 


pA'k p1 

/ d'lpr dx e 
Jo J-1 


i[(n2-Fn3—ni)'0r 


Pmi {x)Pm2 {x)Pm3 (x) , (17) 


— (2?ni ± 1) ^ ^ hni—n3,m2^n3,n 


mi m 2 m 3 

0 0 0 


where we have used the Wigner 3j symbol. 


mi m2 m3 

h h is 


(18) 


(19) 


Once again we can obtain an unbiased estimate of the polarisation 
by calculating this coupling matrix and then inverting it. Explicitly 
the polarisation will be 

Q = 2^{S2,o), (20) 

U = 2^(52, 0 ). (21) 

Just as in Section [±T] where we had to ensure that we included all 
the Fourier modes of the systematics, here we will have to include 
all of the Legendre polynomials that describe f{x). This will de¬ 
pend on the underlying systematics and the range of 'ipt angles seen 
at each pixel. Unlike in Section [TT] where the term chosen is a di¬ 
rect result of the spin of the systematic required to be removed, here 
there is no physical motivation for the terms to use. We simply re¬ 
quire that enough terms are used such that we obtain a satisfactory 
fit for the combined result of the systematics, f{x). 


3 TEST ON SIMULATIONS 


To test the map making algorithms described in Section]^ we sim¬ 
ulate two common types of experiment: one satellite-like exper¬ 
iment, having an extensive range of orientation angles, and one 
balloon-like experiment where this range is limited. To this end, 
we use the EPIC ( |Bock et ^|2009t and the LSPE ( |Aiola et 
|2012| ) scanning strategies respectively. The hit maps of the two scan 
strategies are shown in Fig.|^ In this figure, we also plot the polar¬ 
isation angle coverage (p 2 ) for each pixel. 


P2 


-j -^hits 

J— [cos^( 2 i/>J)+sin^( 2 i/>J)], 

^ ® i=l 


( 22 ) 


where A^hits is the number of hits that a pixel has received. This 
quantity demonstrates the range of 'ipt angles provided by the scan 
strategy. The range of p 2 goes from 0 to 1 and the lower the value, 
the better the polarisation angle coverage. 

In all the simulations we use a fiducial power spectrum with a 
scalar-to-tensor ratio of 0.1 and lensing B-modes are also present. 


3.1 Extensive 'ipt coverage algorithm with a satellite-like 
experiment 

We simulate noisy TODs with systematic errors. The main source 
of systematics we will be considering are leakage from tempera¬ 
ture fluctuations to polarisation fluctuations. Therefore it may be 
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Figure 1. Left panel: The amplitudes of Legendre polynomials required to describe g{'ip) = cos(2'0 + tt/S), where the range is ±0.25 rads. This range 
is typical of the 'ft range of LSPE. Centre panel: A demonstration of using the first 3 Legendre polynomials to describe the systematic g{'f), the crosses are 
the ’’hits” which were randomly generated with a uniform distribution over the full range. The dashed line shows the reconstruction of the systematic. Right 
panel: The error between the reconstruction and the input systematic. The functional form of the residual looks very similar to that of the next polynomial in 
the series, P 3 {x). 



Figure 2. The hit maps (left) and polarisation coverage (right, p 2 ), see equation of the two sc an strategies used in the simulations. The plots show a full 
sky Mollweide projection of the sphere. Upper panels: the EPIC scan strategy jBock et al.|2009) used to demonstrate the “extensive 'ft' range algorithm. 
The EPIC scan strategy is designed to mitigate systematic errors by having many crossing angles for each pixel. This projection is in Galactic coordinates. 
Lower panels: the LSPE scan strategy (Aiola et al.|2012) used to demonstrate the “limited ft” range algorithm. The LSPE gondola will perform scans of 
constant azimuth, changing the elevation ~daily. There is, therefore, a limited range of crossing angles for each pixel. LSPE will obtain good polarisation 
angle coverage using a stepped HWP. This projection is in Ecliptic coordinates. 


sufficient to only simulate the systematics introduced by imperfec¬ 
tions that couple temperature to polarisation. A TOD element, tj , 
is simply the temperature and polarisation response multiplied by 
the underlying CMB sky and then integrated, 

t, = o/du[i,J(„m„ + Ap) + S?(„)0(u + Ap) 

+B7(u)t/(u + Ap)], (23) 

where X(u) is the sky emission in the Stokes parameter X from 
the direction pointed to by the unit vector u. Bf (u) is the beam 
response in the direction u for the Stokes parameter X when ori¬ 
entated in the position j. G is the gain of the detector and Ap is 


the shift in the temperature beam due to the pointing error. The 
position j describes the orientation of the telescope by the stan¬ 
dard Euler angles and the orientation of the HWP. In order to sim¬ 
ulate this correctly one would need to convolve the sky over this 
4 dimensional space. For a high resolution experiment this would 
be computationally infeasible, especially when one requires many 
CMB realisations. We therefore only simulate beam systematics 
that couple temperature to polarisation due to differential elliptic- 
ity. The second and third term of the RHS of equation ( [2^ can be 
calculated simply from a polarisation map of the sky smoothed with 
the axisymmetric component of the beam. In our current model of 
the HWP in CMB experiments, the orientation of the HWP (fh) 
has no effect on the temperature response. Therefore, we only re- 
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quire the convolution of the beam over the 3 dimensional space 
(^, 0, ^t)- This approximation can be formally written as, 

t, = o/du[i,J(um„ + Ap) + _4a(„)Q(u + Ap) 

+4lf(u)t/(u +Ap)], (24) 

where Af(}x) is the axisymmetric component of the beam re¬ 
sponse. The first term of equation ( [24] ) is calculated by a fast pixel 
space convolution code devel oped in | Wallis et al.| ( |2014| ) based 
on the algorithm described in |Mitra et al.| p011j )~The code pro¬ 
duces the temperature field convolved with the asymmetric beam 
as binned in the 3 dimensional space. In the 0 and 4> space we use a 
HEALPix pixellation ( jGorski et al.|2005| ), and in the space we use 
a linear binning. The convolution code calculates the central values 
of the pixels for this 3 dimensional grid. We use Wside=2048 for 
the HEALPix pixelisation and the space is separated into 80 bins. 
The second and third terms of equation ( [24] ) are calculated using the 
SYNEAST program part of the HEALPix package. Each of these codes 
gives us the central values of the pixelised space. We therefore, use 
linear interpolation to calculate the TOD element for a particular 
pointing. 

We use this set up to simulate the TODs for one detector pair 
for a given scan strategy. For the temperature beam, , we use an 
elliptical Gaussian described by 

(25) 

27rqa^ 

Equation ( [25] ) describes the beam for detector 1. The other detector 
has a similar profile except it is rotated by tv/2 to create a differ¬ 
ential ellipticity between the two detectors. We use cr = 3 arcmin 
corresponding to a FWHM of 7 arcmin and the ellipticity parameter 
q = 1.2. 

We include a differential gain between the detectors by simply 
multiplying one detector’s response by a constant gain factor. We 
also simulate a constant differential pointing by simply including 
an offset in one of the detector pointings in our simulation. 

We use the EPIC scan strategy ( jBock et al.|2009| ) in the fol¬ 
lowing simulations with and without a stepped HWP, to simulate 
one detector pair that suffers from the systematics we consider in 
this paper. See Fig.j^for the hit map of the EPIC scan strategy. We 
step the HWP by tt/S every Ihr. For these satellite simulations we 
do not include a Galactic mask. This map-making algorithm works 
in a very similar way to a binned map and only requires the TOD 
data from one pixel to create an estimate of the Stokes Q and U 
of a pixel. It will therefore work equal well regardless of the sky 
coverage. Here we use the entire sky to make the power spectrum 
analysis simple. 


3.1.1 Simulation 1: No HWP included, no noise 

We simulate a noise-free TOD from a detector pair that suffers from 
a differential gain of the two detectors of 1% and a differential 
pointing of 0.1 arcmin which is 1.5% of the 7 arcmin (FWHM) 
beam. We do not include any differential ellipticity because in this 
simulation we do not have a HWP. Without a HWP the spin-2 sys¬ 
tematic created by differential ellipticity cannot be distinguished 
from the spin-2 polarisation signal. Therefore, the technique we 
present in this paper cannot remove the systematic. If an experiment 
needs to remove this systematic we suggest the methods presented 
in |Ade et ^ ( |2014| ) and jWallis et al.| ( |20T4| l. Fig.j^shows the recov¬ 
ered B-modc power spectrum when a simple binned map is made 


from th is TO D compared to one where the algorithm described in 
Section [2 t] is used. We included the terms *§o,o and a§o,i in equa¬ 
tion ([^ to account for the differential gain and pointing. In Fig. [^ 
we can clearly see that the algorithm has removed the bias on the 
recovered B-mode power spectrum as a result of the temperature 
to polarisation leakage. 


3.1.2 Simulation 2: HWP included, no noise 

We simulate a noise-free TOD from a detector pair that suffers from 
a differential gain of the two detectors of 1%, a differential pointing 
of 0.25 arcmin which is 3.5% of the 7 arcmin beam. The differen¬ 
tial ellipticity is created using the beam described by equation (j^. 
Figurej^shows the recovered B-mode power spectrum when a sim¬ 
ple binne d ma p is made and a map using the algorithm described 
in Section [TT] We included the terms a§o,o, *§o,i and *§ 0,2 in equa¬ 
tion (|^ to account for the differential gain, pointing and ellipticity. 
Figureshows that the bias created from the systematics has been 
removed. We also present the result when the same experiment is 
used but with no HWP present and when a binned map is made. 
This illustrates how a HWP can partially mitigate systematics. It 
also demonstrates that even with this mitigation, further systematic 
removal would be required. 


3.1.3 Simulation 3: HWP included, noise included 


We simulate a noisy version of the TOD used in simulation 2 de¬ 
scribed in 3.1.2 We include noise in the TOD of 1 fiKy/s. This is 
optimistic for a CMB experiment. However it is chosen so that the 
recovered B-mode power is easily detected with one detector pair. 
The algorithm can easily deal with noise as the noise simply prop¬ 
agates through the matrix operation to the map in the same way as 
it does in a binned map-making scheme. The algorithm, however, 
does increase the noise on the recovered Q and U measurements 
and it also increases the covariance of the Q and U estimates. This 
increase depends on the scan strategy: the more crossing angles, 
the lower the increase in noise. The limiting cases of this are (1) an 
ideal scan strategy, where the noise increase is zero and (2) where 
the matrix in equation (9} is singular, in which case the effective 
increase in the noise is infinite. With noise simulations we have 
shown that the noise increase for the EPIC scan and this HWP set 
up is 12%. Fig. [^ shows the recovered B-mode power spectrum 
when a simple binned map is made and a map using the algorithm 
described in Section ITT] and where we have removed the noise bias 
in each case. Again we included the terms *§ 0 , 0 , *§ 0,1 and *§ 0,2 in 
equation to account for the differential gain, pointing and ellip¬ 
ticity. Fig. shows that the bias created from the systematics has 
been removed even in the presence of noise. 


3.1.4 Simulation 4: HWP included, no noise, CMB dipole 
included 

The CMB dipole can in principle leak to polarisation through the 
systematics considered in this paper and this effect can be very 
large. Here, we test if the map-making algorithm can remove such 
a level of leakage. For an experiment where the large scale modes 
are not filtered out at the time stream level, the leakage from the 
CMB dipole could be problematic. Experiments of this type do 
have the benefit of being able to calibrate their detectors using the 
CMB dipole ( jPlanck Collaboration et al.|2015bj ). With this benefit 
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Figure 3. We show the recovered S-mode power spectrum for the satellite-like experiment as an example of the map-making algorithm with a extended range 
of angles, see Section [TT] The left column shows the power spectra and the right the fractional error for the simulations. The red dots show the results 
where a binned map is made and the blue the result when our algorithm is used, the input B-mode power spectra is shown in black. First row: the results from 
simulation 1, see Section [3.1.1| where we consider an experiment without a HWP and include differential gain and pointing systematics in the TOD. As there 
is no HWP we cannot apply our method to remove differential ellipticity so it is not included here. Second row: the results from simulation 2, see Section 
|3.1.2| where we consider an experiment with a HWP and included differential gain, pointing and ellipticity systematics in the TOD. The green dots show the 
recovered power spectra where a binned map is made and the TOD was created without considering a HWP to demonstrate the benefit of a HWP at mitigating 
systematics. Third row: the recovered power spectra for simulation 3, see Section [3.1.3| This is similar to simulation 2 but with noise included in the TOD. 
Fourth row: the recovered power spectra for simulation 4, see Section [T.1.4| This is similar to simulation 2 but with the CMB dipole included in the sky model 
and a lowed level of differetial gain, see Section [TL4] The green dots show the recovered power spectra where a binned map is made and the TOD was created 
without considering a HWP to demonstrate the benefit of a HWP at mitigating systematics. 
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in mind one would expect the differential gain for such an experi¬ 
ment to be lower than for ground-based experiments. To reflect this 
effect we lower the level of differential gain in the simulations for 
this section. We simulate a noise-free TOD from a detector pair that 
suffers from a differential gain of the two detectors of 0.2%. The 
other two systematics were kept the same as in simulation 2 in Sec¬ 
tion Figure shows the recovered B-mode power spectrum 
when a simple binnedmap is made and a map using the algorithm 
described in Section [2T| We included the terms *So,o, *§o,i and *§ 0,2 
in equation 12} to account for the differential gain, pointing and el- 
lipticity. Figure shows that the bias created from the systematics 
has been removed. We also present the result when the same ex¬ 
periment is used but with no HWP present and when a binned map 
is made. This illustrates how a HWP can partially mitigate system¬ 
atics. It also demonstrates that even with this mitigation, further 
systematic removal would be required. 

3.2 Limited coverage algorithm with a balloon-like 
experiment 

We test the limited range form of our map making algorithm on 
a balloon like experiment. We use the LSPE scan strategy ( |Aiola| 
|et al.|201^ — see Fig. for the hit map of the LSPE scan strat¬ 
egy. This is a typical balloon-like scan strategy where the gondola 
rotates rapidly to cover the sky. This provides good sky coverage 
(25% for LSPE). However, the pixels are always scanned in a sim¬ 
ilar direction. Although this can make systematic mitigation prob¬ 
lematic, we show here that this problem can be avoided with a suit¬ 
able map-making algorithm. Instead of using the range of crossing 
angles to accurately characterise the systematic and therefore re¬ 
move it, here we assume that the systematics are slow enough func¬ 
tions of 'ijjt that in the small range of probed by the instrument, 
the combined effect can be described by a few Legendre polynomi¬ 
als. The maximum range of 'ijjt for all the pixels in the LSPE scan 
strategy is ?^0.5 rads. Eigure[^shows that with this small range the 
spin-2 systematic (the fastest changing systematic we consider) can 
be recovered to a fraction of a percent with just the first 3 Legendre 
polynomials. 

In the simulation we use the same beam shape as used in|Wal-| 
|lis etal.lpOUl i — see Eigure 1 of that paper, which has a EWHM of 
1.5°. This beam is a simulation of the beam planned to be on board 
LPSE. We simulate TODs from a detector pair — one dectector has 
the same beam as the other but rotated by 7r/2 to provide a differ¬ 
ential ellipticity. LSPE will achieve the required angle coverage by 
using a stepped HWP. We simulate this by rotating the polarisation 
sensitivity of the beam. We step the HWP by tt/S every hour in the 
simulation. We simulate the LSPE scan strategy for 15 days, each 
day the telescope performs scans of constant elevation. We then 
change the elevation daily. 

Unlike in the satellite-like experiment, simulating the asym¬ 
metric beam for polarisation is feasible. Even though the asymme¬ 
try of the polarised beam will in principle create a bias that our al¬ 
gorithm does not remove, the resulting bias is considerably smaller 
than the temperature leakage. The simulation performed here is a 
good demonstration of this as we simulate the systematic errors 
such as differential pointing completely. However, we only remove 
the resulting temperature to polarisation leakage and we do not ac¬ 
count for the resulting polarisation to polarisation errors. 

As the amount of simulated data is much smaller for this low 
resolution balloon-like experiment, we do not have to make the 
same approximations as we do in the satellite-like experiment. We 
perform a pixel based integration of the beam multiplied by the 



Figure 4. We plot the reconstruction of the F^-mode and B-mode power 
spectrum for the balloon-like experiment, see Section [T^ as an example of 
the map-making algorithm for a scan strategy that has a limited range of 
'ipt- We show the recovered power spectrum averaged over 100 realisations 
for maps created in three ways. The diamonds show the recovered power 
spectra for a binned map, the stars correspond to including the term ;S'o,o 
in equation iH- Finally we show the power spectra including both 50,0 
and Bo,i in equation d The resulting recovered power spectra and error 
bars are shown. The bias in the recovered B-mode has been reduced by 2 
orders or magnitude to less than 5% of the error bars, which as described in 
Section [T^ is for an optimistic LSPE noise. 


CMB sky for each TOD element in the experiment. We simulate 
differential gain by multiplying one detector of the pair by 1.01 to 
create a 1% error. Differential pointing is created by changing the 
pointing position of the second detectors beam by 0.05° which is 
3% of the 1.5° FWHM beam. As described above the differential 
ellipticity is created using the simulated beam shown in figure 1 of 
[Wallis et~X] ( |2014] ). We include noise in the TOD corresponding 
to a noise level in the map of 0.1 /xK per 1.5° beam. This level of 
noise is optimistic for LSPE (which should achieve ~ 7 /xK per 
1.5° beam ( jBersanelli et al.|20T^ ). However, as an example of the 
algorithms ability to deal with noise, this level of noise is more than 
sufficient. 

Eigure shows the recovered E- and B-mode power spec¬ 
trum for the balloon-like experiment. As the LSPE scan strategy 
only covers 25% of the sky we use a simple pseudo-C^ estima¬ 
tor ( [Brown et al.|2005| l to recover the polarised power spectra. We 
show the recovered power spectrum averaged over 100 realisations 
for maps created in three ways. Unlike in the extensive range 
case there is little physical interpretation of the terms removed in 
the limited 'ijjt range case. We are creating different approximations 
to the combined effect of the systematic f{x) from equation ( p^ . 
The first map-making algorithm we use assumes /(x)=0. This is 
equivalent to a simply binned map and the recovered power spectra 
are shown as the diamonds in figure The bias from the binned 
map is obvious. We can improve this by including the term Bo,o 
in equation ( p^ . This makes the approximation /(x)=const and 
reduces the bias by an order of magnitude, as shown in Eig. as 
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the stars. Finally we improve this further by assuming that f(x) 
is a linear function. This is done by including both So,o and So,i 
in equation ( p^ . The resulting recovered power spectra and error 
bars are shown in Fig. The bias in the recovered B-mode has 
been reduced by 2 orders or magnitude to less than 5% of the er¬ 
ror bars, which as described above are error bars for an optimistic 
LSPE noise level. 

Figure clearly demonstrates that our technique can remove 
the bias as a result of temperature to polarisation leakage. However, 
this comes at the cost of a noise penalty in the map. Including the 
So,o only and adopting /(a:)=const incurs a 5% increase in the 
noise power with respect to the binned map. The analysis where 
f(x) is modelled as a linear function and both So,o and So,i terms 
are included in equation |T^ creates an increase in the noise power 
of 12% with respect to a binned map. 


4 IDENTIFYING SYSTEMATICS 

The extensive range map-making algorithm relies on being 
able to remove systematics which have different Fourier modes in 
('tpt, V^r) space to that of the polarisation signal. As discussed in 
section 1^ the map-making algorithm requires us to know that the 
systematics are present in the TOD so that we can choose to in¬ 
clude the correct Fourier modes in equation and create maps 
that are clean of those systematics. We stress that we do not need to 
know the exact nature of the systematic. For example, if the exper¬ 
iment is suffering from differential pointing of the detector pair, at 
no point do we need to know by how much or in what orientation 
the beams are misaligned. Moreover, we do not require a temper¬ 
ature map to remove the signal. This is an improvement over the 
method used in |Ade et~^ ( |2014] ). We only need to know that the 
experiment is suffering from differential pointing and therefore we 
know to include the So^i term in the analysis. One down side of 
using this method to clean systematics is the increase in statistical 
noise. Every term included in the analysis increases the statistical 
noise of the recovered polarisation maps, and also the cross corre¬ 
lation of Stokes Q and U. The level of the increase is dependent 
upon the scan strategy — the more extensive {'ipt , ^r) coverage the 
experiment has, the smaller the increase of noise. We demonstrated 
this with simulations of the EPIC scan strategy (see Section|^. The 
increase in the noise power spectrum, going from a binned map 
to our map-making algorithm accounting for all three systematics, 
was 12%. This would be lower if not all systematics were consid¬ 
ered. With this increase of noise in mind it would be undesirable to 
include terms needlessly, but we also do not want to create a bias 
by neglecting a term if the systematic is present. 

We now turn our attention to a practical process to determine 
whether a potential systematic should be removed. We start by 
making a map of the systematic. This is done in the same way that 
the Q and U maps are made. The inversion of equation will give 
us an estimate for the Sn,m terms we included. By taking the real 
and imaginary parts of this we can make a map of our systematics. 
Eigure shows maps of the systematics recovered from simula¬ 
tion 3 described in Section [Tl. 3 1 We also show the predictions for 
each systematic based on prior knowledge, and finally, the differ¬ 
ence between the prediction and recovered systematic maps. This 
figure demonstrates how we can accurately recover the systematics 
without any prior knowledge of the instrumental imperfections or 
of the underlying temperature field. 

The top row of Eigure shows the map of differential gain. 
As this is a spin-0 term there is only one non-trivial map to show. 


The prediction is created by convolving the temperature field with 
the beam and multiplying it by the size of the differential gain. 
The difference is consistent with noise as a result of the noise in 
the TOD. Rows 2 and 3 show maps of differential pointing. Being 
spin-1, row 2 shows the error if the instrument was oriented with 
^t=0 and row 3 shows the error if the telescope was orientated 
with ^t=7r/2. The systematic error on a TOD sample is simply a 
rotation of this spin-1 vector. The prediction is created by using the 
first differential of the temperature map and multiplying the results 
by half the angular size of the differential pointing. Rows 4 and 5 
are maps of the differential ellipticity. As the differential ellipticity 
is a spin-2 systematic we are only able to make the distinction be¬ 
tween the systematic and the polarisation using the HWR The fields 
depicted in rows 4 and 5 are similar to the Stokes Q and U fields 
respectively. The prediction was created by using the underlying 
temperature field and the beam shape used in the simulation. It was 
shown in |Wallis et al.| ( |201^ that the spin-2 systematic would have 
the form. 


'^im 


(26) 

'^im 

■ / T 

(27) 


where Aa^ and Aa^ are the E and S-mode of the system¬ 
atic, aj^ is the spherical harmonic decomposition of the temper¬ 
ature field and is the second azimuthal mode of the spherical 
harmonic decomposition of the difference of the two temperature 
beams in the detector pair. 

Eigure shows that the map-making algorithm can correctly 
recover maps of the systematic errors. However, the maps are noisy. 
One can imagine a situation where the noise level is too large to see 
the systematic but there is a non-negligible effect on the recovered 
B-mode power spectrum. This will be especially true when many 
detectors are considered as the maps would have to be made for 
each detector pair. However, since each of the systematics couple 
to the temperature field, the recovered maps of the systematics will 
correlate with the temperature field. 

In Eig.j^we plot the cross-power spectrum between system¬ 
atic maps and the temperature map created using the TOD used in 
the simulation, where we have also deconvolved for the beam. We 
have over plotted NqCJ, NiiCf, N 2 i‘^Cj', where Ni are normal¬ 
isation factors calculated by minimising the absolute residuals be¬ 
tween the model and the cross-power spectrum. This demonstrates 
the known result that the systematics are simply the temperature 
field, or some derivative of the temperature field. Even in a noisy 
systematic map the cross-power spectrum will provide a valuable 
insight to the size of a systematic effect. We recover a non-zero 
correlation because the systematics are present in the TOD. If this 
were not the case and the TOD was clean then we would find a 
cross-power spectrum consistent with zero. This provides us with a 
recipe to test for the presence or absence of systematic effects: if, in 
an experiment, this cross-power spectrum is shown to be consistent 
with zero then the maps can be (re-)made not accounting for this 
systematic. The resulting increase of noise in the map, with respect 
to a binned map, would thus be kept to a minimum. 


5 DISCUSSION 

We have developed two map-making algorithms to remove com¬ 
mon systematics that couple temperature to polarisation in dif¬ 
ferencing CMB polarisation experiments. The result of the map- 
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Figure 5. Each row shows the systematic from the noisy TOD in simulation 3 described in Section [3.1. 3 1 Each section of the sky is a gnomonic projection 
of a 12.5° square patch of sky. We plot the map extracted using the systematic map making method described in Sectionj^on the left. The middle column 
displays a prediction for the systematic based on prior knowledge. The right column shows the difference which in each case is compatible with noise. Row (a) 
shows the results for differential gain (*So,o), rows (b) and (c), differential pointing (3f^[S'o,i] and ^[*So,i]) and rows (d) and (e) show the differential ellipticity 
systematic ($R[;So, 2 ] and ^[;So, 2 ])- 
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Figure 6. We plot the cross-correlation of the systematic maps (see Fig¬ 
ure with the temperature map created using the TOD in simulation 3 in 
Section]^ The cross-power spectrum is a good method to find systematics 
that may be below the noise level in the systematic maps. The differen¬ 
tial gain systematic error map cross correlated with the temperature map is 
plotted in dark blue dots. The light blue dots show the differential pointing 
maps cross correlated with the temperature map. The differential ellipticity 
map cross correlated with the temperature map is shown as the red dots. 
Over plotted are NqCJ, NilCj, N 2 i^Cj, where Ni are normalisation 
factors. This demonstrates that the systematics are simply the temperature 
field convolved with the beam, or derivatives thereof. 

making algorithms is the polarisation sky smoothed with the ax- 
isymmetric part of the beam used. The systematics we consider are 
differential gain, pointing and ellipticity of two detectors in a de¬ 
tector pair, all of which were shown to be an issue in the BICEP2 
experiment { |Ade et aL|2014| l. The main issue with these systemat¬ 
ics is the leakage from temperature to polarisation that they create. 
[Shimon et al.| ( |2008j ) showed that the coupling from temperature 
to polarisation of these systematics have spin-0,1 and 2 properties 
respectively. We used this understanding to develop the algorithm 
used here. 

The first algorithm, described in SectionjTT] removes the sys¬ 
tematics by separating the Fourier modes of the systematics and the 
Fourier mode of the polarisation, using equation ([^. The technique 
requires a suitable scan strategy. The angle coverage of the orienta¬ 
tion of the telescope, t/;t, must be extensive to allow for the different 
Fourier modes to be distinguished. We have shown through simu¬ 
lations in Section [TT] that the FPIC ( [Bock et al.|2009| ) scan strategy 
provides the required amount of angle coverage. 

In Section [STT] we demonstrated the effectiveness of the algo¬ 
rithm through three simulations. Simulation 1 showed the ability 
of the technique to remove differential gain and pointing when a 
HWP is not used. Without a HWP the spin-2 systematic and the po¬ 
larisation signal are degenerate and, therefore, cannot be separated 
using this technique. In this case we suggest using the methods pro¬ 
posed in [Wallis et al.| ( |201^ and [Ade et al.[ ( [2014] ). The technique 
however, can remove differential gain and pointing as they have a 
different spin to the polarisation signal. In Section [?.1.2[ and [3.1. 3 [ 
we presented simulations including a HWP. In these simulations 
we showed that the technique can simultaneously remove differen¬ 
tial gain, pointing and ellipticity. The CMB dipole can contribute 
to the leakage from temperature to polarisation. We demonstrated 
in Section !?. 1.4[ that the map-making algorithm can deal with this 
level of leakage. 

One draw back to this technique is an increase in the statis¬ 
tical noise in the resultant (cleaned) Q and U maps. The level of 
the noise increase is dependent on the scan strategy — the more 


even the angle coverage, the lower the increase in noise. An 
ideal experiment would suffer no increase in noise. At the other ex¬ 
treme where the coverage is not large enough, the matrix in equa¬ 
tion becomes singular and the effective increase in noise is infi¬ 
nite. Through simulations we have shown that the increase in noise 
power for simulation 3 is 12% when compared to a binned map. 

The second algorithm, described in Section [2?] removes the 
temperature to polarisation leakage by creating a model for total 
leakage as a function of the orientation of the telescope. The com¬ 
bined effect of the systematic is modelled as a smooth function of 
the orientation angle. With this assumption we can then describe the 
combined systematic f{x) by a few Legendre polynomials. Figure 
shows that a spin-2 systematic can accurately be reconstructed 
by the first three Legendre polynomials, with a range of 0.5 
rads. This ^pt range was chosen to be representative of the LSPE 
scan strategy ( [Aiola et al.|2012[ ), where the maximum range is 
?^0.5 rads. In Section [3.2[ we demonstrated that the algorithm can 
remove the temperature to polarisation leakage from differential 
gain, pointing and ellipticity in a simulation where the LSPE scan 
strategy was used with a stepped HWP. As with the extensive 'ipt 
range technique there is an increase in the statistical noise of the 
polarisation maps when using this technique. Through simulations 
we showed that the increase of noise using the LSPE scan strategy 
was 12%. 

In Section [^ we have presented a method to identify if sys¬ 
tematics are present in a TOD. The extensive 'ipt range algorithm 
separates the systematics into Fourier modes and generates an esti¬ 
mate of the polarisation free of these systematics. It also at the same 
time creates an estimate of the systematic. We showed in Fig.[^the 
maps of the systematics recovered from simulation 3 (see Section 
[3.1.3[ ). These maps can be used to identify if a systematic is present. 
However, the noise in the TOD and the relative size of the system¬ 
atic effect could render the reconstructed maps too noisy to see the 
systematic easily. To increase the signal-to-noise, we suggest cal¬ 
culating the cross-power spectrum of these systematic maps with 
the temperature map. As the systematics we are considering are 
due to temperature to polarisation leakage then, if the systematic is 
present, we expect to see a non-zero cross-correlation. We showed 
the cross-power spectrum between the systematic maps and the 
temperature map created using the TOD in Fig.|^ This technique 
can be used to identify if systematics are present. This is crucial 
to test the validity of an experiment’s polarisation maps, but could 
also be used to identify whether a systematic must be removed us¬ 
ing this technique. Accounting only for those systematics that are 
actually present in the TOD would minimise the increase in noise 
associated with the correction algorithms developed in this work. 
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